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Preface 


The  following  report  is  a  revised  version  of  an  invited  peper  "The 
Mathematical  Foundations  of  the  Hotine-Marussi  Approach  to  Geodesy"  which  was 
presented  on  June  5,  1989  at  the  Second  Hotine-Marussi  Symposium  in  Pisa.  The 
talk  in  Pisa  essentially  consisted  of  an  informal  discussion  of  the  basic 
ideas/results  in  the  manuscript  and  was  illustrated  by  examples  and  figures 
not  included  in  the  original  paper.  In  the  present  report  these  examples  and 
figures  are  included  together  with  additional  material  derived  from 
discussions  with  participants  at  the  Symposium.  This  has  resulted  in  almost  a 
50%  increase  in  the  size  of  the  original  manuscript. 

The  material  addresses  the  topological  foundations  which  underlie  the 
entire  Marussi-Hotine  formulation  of  differential  geodesy.  It  also  furnishes 
an  outline  of  some  of  the  mathematical  preliminaries  required  in  the 
investigation  of  the  hierarchy  of  questions  related  to  the  Marussi  Hypothesis 
which  is  the  primary  topic  of  our  research  contract.  A  subsequent  report,  now 
in  preparation,  will  be  devoted  to  a  mathematical  appreciation  of  Marussi 's 
contributions  to  geodesy  and  will  contain  an  introduction  to  the  Marussi 
Hypothesis  in  its  original  formulation. 


1.  Introduction 


The  Scottish  mathematician  and  natural  philosopher  Sir  James  Ivory  began 
one  of  his  papers  in  The  Phi l osop/ti col  Magazine  nod  Journal  for  1824  with  the 
words: 

"It  is  not  my  intention  to  trace  minutely  the  various  labours  of 
philosophers  on  the  Figure  of  the  Earth,  but  to  state  concisely  the 
present  mathematical  theory  on  the  subject,  and  to  add  some 
observations  upon  it." 

Although  I  cannot  aspire  to  attain  either  Ivory's  eloquence,  or  lucidity 
of  exposition,  in  the  present  report  I  will  attempt  to  imitate  his  endeavors 
by  an  examination  of  mathematical  foundations  of  the  Marussi-Hotine  approach 
to  geodesy.  This  approach  is  based  on  an  essential  and  extensive  use  of 
differential-geometric  and  tensor-theoretic  methods  which  abandon  the 

traditional  views  of  H.  Bruns  (1878)  and  F.  Helmert  (1880)  of  geodesy  as  "the 
science  of  the  measurement  and  mapping  of  the  Earth's  surface." 

Indeed,  in  1952  MARUSSI  [1]  summarized  his  approach  by  the  statement  that 

"Geodesy  is  the  science  which  is  devoted  to  the  study  of  the 

Earth's  gravity  field." 

From  this  viewpoint  it  follows  that  our  primary  object  of  concern  is  the 
family  2  of  external  equipotential  surfaces  of  the  Earth,  and  following  F. 
Bocchio  (1972)  we  will  call  the  investigation  of  these  surfaces  differential 
geodesy.  This  terminology  is  slightly  more  general  than  Marussi's  intrinsic 
geodesy  [2],  and  less  general  than  Hotine's  mathematical  geodesy  [3]. 
Although  Marussi's  approach  stressed  the  importance  of  conceptually 

formulating  his  theory  in  a  coordinate- free  manner,  in  practice  much  of  his 
work  was  coordinate-dependent  .  He  assumed  the  existence  of  an  adequate  supply 
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of  natural  physically  relevant  coordinates  which  he  called  intrinsic 
coordinates ,  and  we  term  this  assumption  the  Marussi  Hypothesis .  Hotine  was 
less  insistent  and  emphatic  in  this  regard,  but  much  of  his  analysis  was 
inextricably  tied  to  coordinates.  However,  it  should  be  noted  that  some  of 
his  treatise  foreshadows  the  leg  formalism  of  E.  Grafarend  (1986),  and  B. 
Chovitz  (1982)  has  pointed  out  that  in  his  final  work  Hotine  recognized  the 
efficacy  of  differential  forms.  These  were  subsequently  employed  in  geodesy 
by  Grafarend  (1971,  1975)  and  N.  Grossman  (1974)  and  their  study  of  holonomic 
measurables  has  shed  doubt  on  the  general  validity  of  the  Marussi  Hypothesis. 

Despite  its  visionary  nature  the  approach  of  Marussi  and  Hotine  was 
essentially  of  a  formal  character,  viz  they  did  not  make  any  attempt  to  seek 
either  the  precise  conditions  under  which  their  calculations  were  valid  or  to 
address  questions  of  existence  and  uniqueness.  This  is  not  intended  to  be  a 
criticism,  since  virtually  oil  pioneering  work  which  applies  mathematics  to 
physical  problems  is  formal.  By  necessity,  pioneers  are  more  concerned  with 
the  value  and  utility  of  an  idea,  and  not  with  the  ultimate  consequences  or 
restrictions  under  which  their  methods  retain  their  meaning. 

Our  study  is  focussed  on  the  members  of  the  family  2  and  what  must  be 
assumed  —  both  mathematically  and  physically  —  in  order  for  2  to  be 
susceptible  of  geodetic  and  geometric  study.  In  doing  this  we  make  no  direct 
appeal  to  the  physical  equations  whose  solutions  define  2,  and  consequently 
our  inquiry  may  be  regarded  as  a  geometric  preliminary,  or  a  geometric 
constraint  ,  on  the  theory  of  the  geodetic  boundary  value  problem.  The 
analytical  aspects  of  this  problem  have  been  extensively  investigated  by 
Bjerhammar,  Holota,  Hormander,  Moritz,  Sacerdote,  Sansd,  and  Tscherning. 

Our  immediate  concern  is  with  the  general  properties  of  surfaces  in  2, 
and  these  include  not  only  local  properties  where  one  deals  with  a  piece  of  a 
equipotential  surface  S,  but  glolxil  properties  characterizing  the  entire 
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surface  S.  We  will  see  that  both  of  these  require  topological  notions. 
Indeed  if  differential  geometry  is  regarded  as  providing  a  geometric 
constraint  on  questions  in  physical  geodesy,  then  topology  furnishes  a 
topological  constraint  on  differential  geometry. 

Topological  considerations  are  perhaps  new  to  mathematical  geodesists; 
however,  they  have  played  a  significant  role  in  modern  mathematics.  They  have 
allowed  mathematicians  to  systematically  dissect  and  analyze  the  fundamental 
notions  of  what  is  meant  by  the  terms  space,  continuity,  differentiability, 
etc.  This  has  permitted  one  to  recognize  the  aspects  of  these  concepts  which 
are  strictly  Euclidean  as  well  as  those  which  are  more  general.  Indeed,  in 
1939  Hermann  Weyl,  one  of  the  giants  of  twentieth  century  mathematics, 
remarked  that 

"In  these  days  the  angel  of  topology  and  the  devil  of  abstract 
algebra  fight  for  the  soul  of  each  individual  mathematical  domain." 
The  last  half-century  of  mathematical  endeavors  have  vigorously  confirmed 
Weyl's  comments. 

Mathematically  such  considerations  are  not  new,  and  date  back  to  Leibniz 
(1678),  who  proposed  their  investigation  as  analysis  situs.  This  name  was 
employed  into  the  early  part  of  the  twentieth  century,  when  the  term  topology 
came  into  popular  use.  It  is  amusing  to  note  that  the  same  man  who  coined  the 
term  geoid  in  1873  invented  the  word  topology  in  1847.  He  also  published  an 
account  of  the  one-sided  Mobius  strip  four  years  before  A.F.  Mobius  (1865)1 
This  was  Gauss'  student  J.B.  Listing  (1808-1882),  and  his  "Vorstudien  zur 
Topologie"  was  one  of  the  earliest  systematical  treatments  of  the  subject. 
Listing's  view  of  topology  as  a 

"kalkulatorische  Bearbeitung  der  modalen  Seite  der  Geometrie," 
is  in  a  sense  how  we  will  seek  to  apply  it  to  differential  geodesy. 

In  differential  geodesy  tu>o  important  types  of  assumptions  arise  when  one 
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attenpts  to  pass  beyond  the  formal  aspects  of  the  theory.  We  will  call  them 
the  vicinity  and  smoothness  assumptions  and  show  that  they  embody  mathematical 
and  physical  requirements. 


2.  The  Vicinity  Assumption 


The  uicinity  assumption  incorporates  at  least  two  ideas.  Physically,  it 
begins  with  the  idealization  of  the  figure  of  a  uniformly  rotating  Earth  being 
roughly  spherical.  More  precisely,  one  regards  the  Earth  as  being 
topologically  equivalent  to  a  sphere  which  includes  the  shape  being  that  of  a 
sphere,  a  spheroid,  or  an  ellipsoid.  One  then  asks  under  what  circumstances 
can  one  take  the  surfaces  of  I  to  be  closed  surfaces?  Intuitively,  the 
answer  is  obvious:  the  equipotential  surfaces  must  lie  'near'  the  surface  of 
the  Earth.  However,  a  much  more  precise  answer  has  long  been  known  —  but 
almost  forgotten! 

The  answer  was  given  in  1901  by  Paolo  Pizzetti  (1860-1918),  who  from 

1900-1918  was  Professor  of  Higher  Geodesy  at  the  University  of  Pisa.  In 

PIZZETTI  [4],  he  proved  that  when  the  maximum  and  mean  radii  R  ,  R  of  the 
r  max  m 

Earth  are  subject  to  the  inequality 


(11 


R  -  R  <  R  , 
max  m  100  m 


with  ratio  m  of  the  centrifugal  and  gravitational  accelerations  being 

1 


(2) 


m  < 


280  ' 


then  the  surfaces  of  1  are  closed  for  altitudes  h  above  the  surface  of  the 

Earth  which  satisfy  the  Pizzetti  inequality: 

(3)  0  <  h  <  5R  . 

m 

He  also  observed  that  (3)  remains  valid  when  lunar  tidal  effects  are  involved, 
i.e.  m  <  The  coefficient  of  the  upper  bound  in  (3)  to  4-place  accuracy 

is  4.8930,  but  of  course  Pizzetti 's  values  for  the  geophysical  parameters  are 
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now  obsolete.  However,  upon  taking  latest  accepted  I.A.G.  values  given  by 
Moritz,  (1988),  with  m  <  to  3-place  accuracy  this  coefficient  is  4.969. 
The  choice  of  the  coefficient  5  in  (3)  is  of  no  particular  significance 
except  it  works!  The  essential  idea  is  that  on  the  basis  of  very  general 
assumptions  PIZZETTI  was  able  to  establish  such  an  inequality. 

In  (5)  VJ.D.  Lambert  presented  an  excellent  qualitative  discussion  of  a 
similar  situation,  and  traced  such  considerations  back  to  P.S.  Laplace  and  his 
"Traite  de  M^canique  Celeste,  Livre  III"  (see  Chapter  VII  of  the  Bowditch 
translation  of  1832).  Laplace  analyzed  the  case  of  the  equipotential  surface 
'marking  the  outermost  limit  of  the  atmosphere'  and  determined  its  equatorial 
radius  to  be  6.6  times  the  radius  of  the  Earth.  Lambert  did  not  make 
reference  to  the  work  of  Pizzetti. 
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Figure  1.  The  gravitational  potential  V  is  illustrated  in  the  left  figure, 
the  gravity  potential  W  in  the  right  figure.  The  corresponding 
equipotential  surfaces  for  V  and  W  are  indicated  by  solid  and  dashed 
lines.  This  illustration  is  taken  from  page  537  of  an  article  by  h.  JUNG  in 
the  Handbuch  der  Physik,  Bd.  47  (Springer-Verlag,  1956). 
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It  is  immediate  that  when  (3)  holds  the  family  2  is  not  only  closed, 
but  also  bounded.  If  we  restrict  our  considerations  to  non-relativistic 
effects,  and  assume  the  space  in  the  neighborhood  ot  the  Earth  is 
3-dimensional  and  Euclidean,  then  each  S  <r  2  is  a  compact  surface.  We  will 
return  to  this  important  topological  property  in  a  moment.  Henceforth  we  will 
consider  only  such  S  €  2. 

The  second  aspect  of  the  vicinity  assumption  is  to  emphasize  that  almost 
all  of  the  analysis  involving  differential  geometric  and  tensorial  techniques 
is  local  and  valid  only  on  a  piece,  viz  in  the  vicinity  of  a  point  P,  of  a 
surface  S  €  2.  Stated  more  precisely,  these  considerations  hold  only  on  an 
open  connected  subset,  i.e.,  a  domain,  of  S.  When  this  is  provided  with  a 

r 

coordinate  system  x  this  domain  becomes  a  coordinate  neighborhood ,  or  chart 
about  the  point  P.  Hence,  in  general,  coordinates  are  local  in  character, 
i.e.,  valid  only  in  a  domain  of  S,  and  in  order  to  obtain  a 
"coordinatization"  of  even  a  piece  of  S  it  is  necessary  to  employ  several 
charts.  The  charts  must  be  compatible,  i.e.  "mesh-together,"  and  we  will 
return  to  this  in  Section  3.  These  charts  are  said  to  cover  the  piece  of  S, 
and  by  virtue  of  the  compactness  of  each  S  e  2  the  Heine-Borel-Lebesgue 
covering  property  asserts  that  the  entire  surface  S  can  be  covered  by  a 
finite  number  of  charts.  This  is  illustrated  by  taking  S  to  be 
2-dimensional  sphere  S£  where  two  charts  are  required  to  "patch  up"  the 
coordinate  singularities  which  occur  at  the  North  and  South  poles.  An  even 
better  covering  of  S2,  which  avoids  equatorial  singularities,  employs  six 
charts.  In  Section  4  we  will  see  that  such  singularities  are  not  the  result 
of  an  inappropriate  choice  of  coordinates,  but  a  topological  property  of  S2 
and  surfaces  which  are  topologically  equivalent  to  it. 
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Figure  2.  This  figure  is  entitled  "Limiting  surfaces  of  the  atmosphere  and 
adjacent  level  surfaces  ..."  and  is  taken  from  LAMBERT  [5).  Solid  lines 
denote  limiting  equipotential  surfaces,  dashed  lines  denote  analytical 
continuations  (for  a  detailed  description  see  [5]). 


3.  The  Smoothness  Assumption 


The  smoothness  assumption  refers  to  the  degree  of  differentiability  of 
the  functions  and  coordinates  occurring  in  our  analysis.  We  briefly  review 
the  relevant  notions  for  functions  defined  on  a  domain  ft  of  a  3-dimensional 
Euclidean  space  E^. 

A  function  F  defined  on  ft  is  said  to  be  of  class  C  for  an  integer 
k  >  0  whenever  F  and  all  its  partial  derivatives  up  to  and  including  those 

of  order  k  are  continuously  differentiable  with  respect  to  the  coordinates 

r  k  k 

x  on  ft.  In  this  case  we  write  F  e  C  (ft),  or  simply  F  €  c  when  ft  is 

evident.  The  class  corresponds  to  merely  continuous  functions  and  &' 

indicates  (real)  analytic  functions.  In  general  0  <  k  <  ®  <  w  and  if  no 

specific  value  k  >  0  is  indicated,  it  is  convenient  to  merely  say  that  the 

functions  are  smooth.  An  important  generalization  of  these  ideas  was 

introduced  by  E.  Holder  (1882)  in  his  classical  work  on  potential  theory,  and 

we  will  also  employ  the  Holder  class  Ck,?  for  k  >  0  and  0  <  P  <,  1  in 

which  the  kth  order  partial  derivatives  satisfy  a  uniform  Holder  condition  of 

k  P 

order  P.  If  0  <  P  <  1  membership  in  C  '  is  weaker  than 
differentiability,  but  stronger  than  mere  continuity. 

Quite  obviously  the  first  aspect  of  the  smoothness  assumption  is  to 
insure  that  all  the  functions  belong  to  the  appropriate  classes  in  order  that 
the  various  operations  performed  on  them  are  meaningful.  These  requirements 
are  not  stated  in  (1),  [2]  or  [3],  and  often  they  are  likewise  taken  for 
granted  in  texts  on  differential  geometry  and  tensors.  Such  a  naive  and 
uncritical  approach  is  often  harmless,  but  in  the  case  of  surface  theory  one 
clearly  wants  S  to  admit  continuously  varying  tangent  planes,  normal 
vectors,  principal  curvatures,  etc.  and  none  of  these  properties  is  trivial  or 
automatically  true.  These  considerations  will  be  discussed  further  in  Section 
4. 
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An  even  more  important  aspect  of  the  smoothness  assumption  is  to 
ascertain  when  the  data  obtained  from  geodetic  measurements  is  sufficiently 
smooth  to  allow  us  to  even  meaningfully  speak  of  a  piece  of  S.  The 
irregularities  in  the  equipotential  surfaces  of  I  are  dramatically  shown  in 
the  beautiful  computer-generated  pictures  recently  given  by  P.J.  MELVIN,  [6]. 
While  we  may  assume  our  current  data,  or  our  extrapolations/approximations  of 
it  for  a  S  €  2,  are  adequate  for  an  application  of  differential  geometric 
techniques,  as  a  first  step  we  should  unequivocally  establish  the  precise 
nature  of  these  requirements. 

4.  Local  Surface  Theory 

The  notion  of  a  smooth  S  €  1  in  is  very  complicated.  Technically 

it  is  that  of  a  2-dimensional  smooth  oriented  Riemannian  manifold 

isometrically  imbedded  in  .  The  process  of  developing  the  requisite  ideas 

is  lengthy.  It  begins  with  a  connected  topological  space  X  (consisting  of 

suitable  collections  of  sets  called  a  topology)  in  which  the  elements  of  the 

sets  are  point s.  One  next  inposes  axioms  of  separation  and  separability  (the 

Hausdorff  and  Second  Countability  Axioms)  in  order  to  obtain  abstract 

Euclideanlike  properties.  Coordinates  are  introduced  by  making  X  into  a 

topological  manifold,  and  the  dimension  of  the  manifold  is  obtained  as  a 

byproduct  of  this  procedure.  Locally,  a  topological  manifold  Xn  looks  like 

a  domain  of  E^.  The  apparatus  of  differential  calculus  is  obtained  only  when 

X^  is  provided  with  a  differentiable  structure  and  becomes  a  differentiable 

manifold  Mn«  Finally,  Mn  must  be  supplied  with  a  smooth  Riemannian 

structure ,  given  an  orientation,  and  the  resulting  smooth  Riemannian  manifold 

V  imbedded  in  an  Euclidean  E  . 
n  m 

In  his  report  to  the  Third  Hotine  Symposium  (Siena,  1975),  GROSSMAN  (7) 
gave  a  delightful  discussion  of  the  possibilities  and  difficulties.  One  of 
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his  comments  relative  to  the  above  procedure  bears  repeating  since  it  is  just 
as  true  now  as  it  was  then.  He  said 

"The  predominant  view  is  that  space  near  the  Earth  is  a  manifold. 
The  unpleasant  fact  of  life  is  that  no  one  has  ever  described  a 
method  for  coordinatizing  that  supposed  manifold  in  a  way  consonant 
with  physical  reality." 

Stated  more  succinctly  we  can  paraphrase  this  by  saying  that  although  the 

procedure,  or  predominant  view,  is  clear  it  remains  an  open  challenge  to 

implement  it  in  practice.  Until  this  is  done,  much  of  differential  geodesy  is 

at  best  wishful  thinking,  and  at  worst  a  mathematical  fiction. 

Keeping  Grossman's  comment  in  mind,  let  us  survey  the  situation  in  view 

of  the  vicinity  and  smoothness  assumptions. 

The  first  step  is  to  introduce  a  Gaussian  parametrization  of  the 

Euclidean  coordinates  xr  on  a  piece  of  S.  By  the  Implicit  Function  Theorem 

12  3 

this  amounts  to  solving  the  implicit  equation  of  S  :  F(x  ,  x  ,  x  )  =  0  for 
one  of  the  coordinates  and  introducing  an  arbitrary  set  of  parameters 

CT  1  2 

u  =  (u  ,  u  )  into  the  coordinate  chart  on  S.  For  those  S  which  are 
topologically  equivalent  to  a  2-sphere  S2  this  entails  a  square  root 
involving  the  other  two  coordinates,  and  choosing  the  two  signs  on  this 
radical  yields  two  coordinate  charts  on  S.  The  above  process  can  be  done  in 
three  different  ways  (one  for  each  x  )  and  this  leads  to  the  six  charts 
mentioned  previously. 
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Figure  3.  The  covering  of  the  2-sphere  S2  by  six  charts  as  given  in  the 

textbook  Differential  Geometry  of  Curves  and  Surfaces  by  M.  DO  CARMO,  page  56 
( Prentice-Hall ,  1976) .  For  more  details,  lee  our  Appendix,  equations 
(A.l)— (A.3). 


Suppose  this  has  been  done,  and  that  a  piece  of  S  is  locally  given  a 

Gaussian  parametrization 

i  a\  r  r,  1  2. 

(4)  x  -  x  (u  ,  u  )  . 

Then  the  u“  are  regarded  as  independent  parameters  on  the  coordinate  chart 

r 

(U,  x  )  of  S,  and  this  is  said  to  be  a  C  (k  >  0)  piece  of  S  whenever 

IT  k  ,v 

the  coordinates  x  in  (4)  are  C  functions  of  the  parameters  u  .  This  is 

•  r  k  j( 

written,  x  €  (T  and  S  is  locally  a  C  -surface  whenever  the  coordinate 

charts  covering  S  are  selected  in  a  manner  such  that  a  common  value  of  k 

can  be  assigned  to  all  the  charts  covering  S.  As  mentioned  in  Section  2,  by 

the  compactness  of  each  S  €  I,  a  finite  number  of  charts  suffices  to  cover 

each  S  €  2.  Let  ua  and  u1  denote  the  parameters  in  a  pair  of  charts 
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whose  intersection  U  fl  U  contains  an  arbitrary  point  P  €  S. 
and  u1  are  related  by  the  parameter  changes 

—a  .  1  2 . 

u  =  f  (u  ,  u  )  , 


Then  the  u 


(5) 


a  n  .  —  1  —2  . 

U  =  g  (u  ,  u  )  , 


where  f”,  gn  €  c^.  In  other  words,  the  local  parameter  descriptions  are 
reversible  and  C  -related  to  each  other.  in  order  to  insure  that  outward 
pointing  unit  normal  vectors  n  on  S  retain  their  outward  pointing 


character,  it  is  necessary  that  S  be  orientable.  This  requires  that  the 
functional  determinants  have  a  fixed  sign,  viz  an  orientation,  say 


(6) 


1  2, 


>  0  , 


3(u  ,u  ) 

and  this  guarantees  reversibility  of  the  parametrizations  in  (5). 

K  r 

The  C  -character  of  (5)  is  inherited  by  the  coordinates  x  in  (4).  The 

Inverse  Function  Theorem  then  requires  that  of  least  one  of  the  functional 

determinants 


m  3(x1,x2)  d(x2,x3)  cKxW) 

d(u  ,u  )  cJ(u  ,u  )  5(u  ,u  ) 

be  non-vanishing,  and  this  is  tantamount  to  the  linear  independence  of  the 

r  1  r  2 

tangent  vectors  dx  /du  ,  dx  /du  on  the  piece  of  S,  i.e.,  to  the 
non-vanishing  of  a  component  i>r  of  n.  If  all  three  of  the  determinants  in 

|f 

(7)  vanish  then  the  Inverse  Function  Theorem  fails  to  hold,  and  the  x  are 
said  to  have  a  singularity.  Some  of  these  may  be  artificial,  i.e.,  due  to  a 
poor  choice  of  parameters,  while  others  are  of  a  topological  nature  and  cannot 
be  avoided.  This  is  responsible  for  the  non-existence  of  a  single  global 
singularity-free  coordinate  system  on  S2,  and  explains  the  necessity  of 
several  local  coordinate  charts  on  S2  and  its  topological  equivalents.  In 
effect,  ultimately  this  is  a  byproduct  of  a  celebrated  theorem  of  L.E.J. 
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Brouwer  (1909-10),  see  J.  MILNOR,  (8)  and  (9),  that  an  n-sphere  Sn  admits  a 
globally  defined  field  of  non-vanishing  smooth  tangent  vectors  if  and  only  if 
n  is  odd.  We  will  return  to  this  question  at  the  end  of  Section  5. 

The  above  discussion  has  indicated  the  meaning  and  delicacy  attached  to 
the  requirement  that  S  be  a  C  -surface.  In  the  remainder  of  this  section  it 
is  understood  that  references  to  S  are  strictly  local  in  character,  and 
refer  only  to  a  piece  of  S.  We  now  determine  the  smoothness  requirements  on 
the  various  geometric  quantities  associated  to  S  in  Gaussian  differential 
geometry. 

Mathematically,  the  most  striking  property  of  S  is  that  it  is  imbedded 
in  E^ .  Let  xr  be  Cartesian  coordinates  on  a  domain  Q  c  E^,  with  6^ 
being  the  components  of  the  constant  Euclidean  metric  tensor,  i.e.,  the 
"infinitesimal"  square  of  the  distance  between  a  pair  of  points  P,  Q  €  (7  is 
given  by 

(8)  ds2  -  6  dxrdxs  . 

t  5 

If  now  ua  are  local  parameters  on  a  surface  S  in  E^  with  P,  Q  e  S,  then 
we  also  have  (4)  and  for  P,  Q(  (in  S: 


(9)  ds2  =  a  nd\T dvP 

afi 

where  a ^  are  the  components  of  the  surface  metric  tensor,  i.e.  the 
coefficients  of  the  first  fundamental  form .  The  agreement  of  (8)  and  (9)  on 
n  n  s  immediately  yields 


(10) 


dxr  3xS 


oft 


rs 


du 


r  ex P 


which  is  not  only  the  defining  equation  for  the  a^,  but  also  the  Ixisic 

2 

imbedding  equation.  The  equality  of  these  two  expressions  for  ds  is  the 
identity  of  the  notions  of  distance  on  E^  and  S.  This  is  what  is  meant  by 
an  i some  trie  imbedding  of  S  in  E^ .  It  is  a  nontrivial  requirement  since  it 
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involves  S  being  a  subspace  of  E^  and  possessing  both  a  metric  and 
topology  induced  by  those  of  .  For  the  moment,  let  us  suppose  that  this 
imbedding  has  been  achieved  —  we  will  consider  it  in  more  detail  at  the  end 
of  this  section. 

Since  S  is  a  Ck-surface  (k  >  1)  by  definition  the  xr  e  Ck,  and  hence 

k  1 

by  virtue  of  (10)  it  immediately  follows  that  the  a  e  C  .  Moreover, 

since  the  determinant  a  =  detlla^pll  is  a  polynomial  function  it  follows  that 

a,  Va,  and  the  Levi  -Ci  vi  ta  surface  dualizors  e  are  of  class  Ck 

aP 

By  inspection  of  their  definitions,  the  components  ur  of  the  unit  normal  u 

to  S  are  also  of  class  Ck_1,  and  e  Ck~2 .  The  coefficients  b^  of 

the  second  fundamental  form,  the  determinant  b  =  detllb  „ll ,  as  well  as  the 

ap 

Gaussian  (total)  curvature  K,  the  Germain  (mean)  curvature  H,  and  the  higher 

k-2 

order  fundamental  tensors  (see  [10])  are  all  of  class  C  .  Hence,  when 
k  >  3  all  the  basic  quantities  occurring  in  surface  theory  are  smoothl 

Such  a  smoothness  requirement  allows  one  to  answer  the  question  of  the 
local  existence  and  uniqueness  of  a  surface  S  which  is  assumed  to  be 

isometrically  imbedded  in  E^.  0.  Bonnet  (1867)  proved  the  following: 

Fundamental  Theorem  of  Local  Surface  Theory 

Let  the  a^p,  b^  be  defined  on  a  chart  of  S  c  E^.  Then  a 
necessary  and  sufficient  condition  that  a  pair  of  quadratic 
differential  forms  having  these  tensors  as  coefficients,  be  the 
first  and  second  fundamental  forms  of  S  is  that  locally  S  be  of 
class  ,  with 

(i)  a  >  0  , 
and  that  these  tensors  satisfy 

(ii)  the  Gauss  equation 

R  „  =  b  .  b„  -  b  b„.  , 

a) 3Xp  aX  Pi j  aji  PX 
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and 


( iii )  the  Codazzt -Malnardi  equations 

b  _  -  b  _  =  0  . 

aft, y  cn  ,P 

Moreover,  if  (i)— (iii)  are  satisfied,  then  locally  a  C^-surface  S 
is  uniquely  determined  up  to  position,  i.e.,  a  rotation  and 
translation,  in  E^. 

Both  (ii)  and  (iii)  are  integrability  conditions  of  a  system  of 
differential  equations  known  as  the  Gauss  formulae,  and  (ii)  is  essentially 
the  Theorema  Egregium  of  Gauss.  An  obvious  question  is  whether  Bonnet's 


theorem  can  be  improved,  i.e.,  can  the  differentiability  conditions  be 

2 

weakened?  In  other  words,  can  these  results  be  established  for  C  -surfaces? 


This  would  be  attractive 
ur  €  C1,  and  b  €  C°, 


since  in  such  a  situation:  xr  €  c^,  a 


which  would  still  allow  for  a  continuously 


differentiable  tangent  plane  and  normal  t>.  This  was  done  by  P.  Hartman  and 


A.  Wintner  (1950)  who  by  an  ingenious  argument  replaced  (ii)  and  (iii)  by 
integral  equations.  No  other  improvement  of  the  theorem  is  known,  or  seems 
possible. 

We  now  return  to  the  question  of  the  existence  of  a  local  isometric 
imbedding  of  S  in  E^.  By  definition  this  involves  finding  a  singular!  ti/* 
fiee  solution  xr(u\  u^ )  of  (10),  which  is  said  to  provide  a  realization  of 


a^.  This  requires  regarding  (10)  not  as  a  definition  for  a^,  but  as  a 
system  of  three  non-linear  first  order  partial  differential  equations  for  the 
x  when  the  a^  are  specified.  The  obvious  procedure  is  to  appeal  to  the 
well-known  Cauchy-Kovalevska  Theorem;  however  this  requires  S  to  be  a  (real) 
analytic  surface,  i.e.,  xr  and  a^  €  cf*’.  This  is  a  severe  requirement,  but 
unfortunately  a  common  one  in  pure  mathematics! 


Suppose  for  the  moment  we  generalize  the  situation  and  replace  S  by  an 
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n-dimensional  Riemannian  manifold  V  ,  and  E-,  by  E  .  Then  the  Greek 

n  i  ■*  m 

indices  range  over  the  values  1  to  n,  and  the  Latin  indices  over  1  to  m 
in  (10).  L.  Schlafli  (1873)  conjectured  that  in  this  case,  since  aQ^  has 
n(n  +  1 )/2  independent  components,  it  suffices  to  take  m  =  n(n  +  l)/2.  This 

was  proven  by  M.  Janet  (1926),  and  some  missing  details  were  later  supplied  by 

C.  Burstin  (1931).  Hence  we  have 

The  Fundamental  Theorem  on  Local  Isometric  Imbedding 

Any  Riemannian  manifold  having  a  C^-metric  tensor  can  be 

locally  isometrically  imbedded  analytically  in  an  Em  with 

m  =  n(n  +  l)/2. 

Hence,  when  «  S,  so  that  n  *  2,  the  problem  is  solved  for 

(^-surfaces  by  taking  m  -  3.  It  has  proved  remarkably  difficult  to  replace 

the  stringent  requirement  of  analyticity  by  weaker  conditions  of 

differentiability.  The  best  result  in  this  direction  is  due  to  H.  JACOBOWITZ 

(11],  who  proved  the  following: 

Weakened  Fundamental  Theorem  on  Local  Isometric  Imbeddings 

k  P 

Let  S  be  a  surface  which  admits  a  metric  a  „  e  C  '  with 

aP 

k  P 

k  >  3  and  0  <  C  <  1,  and  which  possesses  a  local  C  ' 

realization  in  E  .  Then 
m 

(i)  if  S  has  Gaussian  curvature  K  >  0,  then  m  =  3,  and 

(ii)  with  no  restrictions  on  the  sign  of  K,  then  m  =  4. 

It  is  likely  that  this  is  the  best  possible  result,  since  A.V.  Pogorelov 

2  1 

(1971)  constructed  a  S  having  a  C  '  metric  a^  which  admits  no 
2 

C  -realization  in  E^.  The  classic  example  is  that  of  D.  Hilbert  (1901),  viz. 
a  (^-surface  having  constant  K  *  -1,  which  cannot  be  analytically  realized  in 
without  the  occurrence  of  singularities. 

The  above  discussion  is  intended  to  illustrate  that  the  question  of  a 
local  isometric  imbedding  of  S  in  is  by  no  means  either  obvious  or  even 
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possible. 


5.  Topological  Considerations 

In  Section  4  we  have  indicated  the  steps  involved  in  the  notion  of  a 
smooth  surface  S  €  2  which  is  imbedded  in  an  E^ .  We  now  discuss  how 
various  global  topological  notions  enter  into  the  picture. 

is  a  topological  space,  and  its  usual  topology  t  is  Hausdorff  and 
possesses  a  countable  basis.  We  denote  this  topological  space  by  (E^,  t). 
Since  S  is  topologically  a  subs/iace  of  ( ,  r),  its  induced  topology  is 
also  Hausdorff  and  second  countable.  By  Pizzetti's  result,  each  S  is  closed 
and  bounded,  hence  since  2  is  in  E^ ,  each  S  e  2  is  com/xict  . 

E^  is  also  a  3-dimensional  topological  manifold,  and  each  S  is  a 

2-dimensional  topological  submanif old  of  E^ .  Moreover,  both  E^  and  each  S 
are  metric  spaces,  which  we  denote  by  ( ,  p)  and  (S,  o ) ,  where  p  and  a 
are  respectively  the  metrics,  i.e.  distance  functions,  on  E3  and  S.  The 
metrics  p  and  a  then  generate  the  usual  topologies  on  E^  and  S,  with  a 
being  induced  by  p.  It  can  be  shown  that  E^  and  each  S  can  be  provided 
with  compatible  differentiable  structures  and  hence  E^  and  S  are 
differentiable  manifolds.  When  provided  with  the  Riemannian  metrics  6 and 
(the  former  being  regarded  as  a  flat  Riemannian  metric),  both  E^  and 
each  S  are  Riemannian  manifolds.  As  a  metric  space,  the  distance  a(P,  Q) 
between  a  pair  of  points  P,  Q  €  S  is  given  by  the  in  ft  mum  of  lengths  L(C) 
of  all  piecewise  smooth  curves  joining  P  and  Q. 

Our  immediate  goal  is  to  show  that  each  (S,  o)  is  a  complete  metric 
space.  By  the  property  of  completeness  we  mean  that  every  Cauchy  sequence  of 
points  on  S  is  a  convergent  sequence.  Stated  intuitively,  completeness 
guarantees  that  every  sequence  of  points  on  S  which  tries  to  converge 
succeeds  in  the  sense  that  it  finds  a  point  in  (S,  a)  to  converge  to.  In  a 
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moment  we  will  translate  this  into  a  geometric/geodetic  property  of  S. 

The  completeness  of  (S,  o)  is  easily  established.  First,  a  standard 
theorem  in  topology  states  that  a  compact  metric  space  is  automatically  a 
complete  metric  space.  An  alternate  argument  is  that  the  metric  space 
( E^ ,  p )  is  complete,  and  that  a  subspace  (S,  n)  of  ( E^ ,  p)  is  complete 
whenever  S  is  closed  in  .  But  this  is  guaranteed  by  Pizzetti's  result! 
Hence  the  metric  space  (S,  a)  is  complete. 

We  now  assume  that  each  (S,  o)  is  a  connected  space.  This  means  that 
each  S  consists  of  a  single  piece,  i.e.  it  is  not  the  union  of  a  pair  of 
disjoint  pieces.  Such  a  requirement  is  obvious  physically,  for  if  it  were  not 
true  there  would  be  no  possibility  of  moving  from  point  to  point  on  S.  The 

geometric/geodetic  consequences  of  completeness  are  then  given  in  the 

following  result  of  H.  HOPF  and  W.  RINOW  [12]. 

Hopf-Rinow  Theorem 

If  a  connected  surface  S  is  complete  as  a  metric  space,  then 

any  two  points  P  and  Q  of  S  can  be  joined  by  a  unique 

geodesic  /  which  has  length  o(P,  Q)  and  is  the  curve  of  shortest 
length  joining  these  points. 

Indeed  a  classical  result  of  0.  Bonnet  (1855)  goes  even  further,  and 
provides  a  bound  on  o(P,  Q): 

Bonnet's  Theorem 

On  a  compact  surface  S  having  Gaussian  curvature  K, 
0  <  kQ  <  K,  where  kg  is  a  constant,  the  distance  is  bounded  by 
rr(P,  Q)  <  rr/Vk^  . 

In  [1]  Marussi  stated  that  a  goal  of  his  intrinsic  geodesy  was  to  carry 
out  geodetic  investigations  without  imposing  any  additional  hypotheses 
relative  to  the  selection  of  a  particular  spheroid  or  ellipsoid.  The  Hopf- 
Rinow  theorem,  which  essentially  involves  translating  a  geometric/geodetic 
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situation  into  a  topological  setting,  assumes  us  that  geodesics  on  the 
equipotentia!  surfaces  S  €  1  are  indeed  curves  of  minimal  length l 

We  are  now  in  a  position  to  obtain  a  rather  sharp  global  topological 
characterization  of  all  S  e  I.  Assume  that  each  S  has  K  >  0  (with  K 
not  necessarily  constant)  or,  equivalently,  that  each  S  is  a  convex  surface. 
Geodetically,  such  surfaces  are  characterized  by  the  fact  that  on  them  plumb 
lines  (i.e.  the  outward  pointing  normals)  at  neighboring  points  do  not 
intersect.  By  cur  previous  analysis  we  know  that  each  S  is  a  complete 
surface.  In  E^  we  then  have  the  following  classification  theorem  of  A.D. 
ALEKSANDROV  [13]. 

Aleksandrov  Theorem 

EVery  complete  convex  surface  S  in  E^  is  topologically 
equivalent  to  either  a  2-sphere  S^,  a  plane,  or  a  cylinder. 

The  latter  two  possibilities  have  K  =  0,  and  hence  can  be  eliminated. 
Thus  all  our  S€  2  are  topologically  equivalent,  i.e.  homeomorphf  c ,  to  a 
2-sphere  S2.  This  equivalence  includes  spheres,  spheroids,  ellipsoids,  as 
well  as  "bumpy"  spheroids,  and  each  of  these  satisfies  the  topological 
properties  of  being  compact  and  connected  as  well  as  being  complete  and 
orientable.  Obviously  S2  is  trivially  a  topological  space,  a  metric  space, 
a  topological  manifold,  and  is  a  subspace  of  E^.  Moreover,  S2  can  be 
provided  with  a  unique  differentiable  structure,  a  smooth  Riemannian  metric, 
and  as  a  Riemannian  manifold  it  is  isome'.rically  imbedded  in  E^.  Hence  each 
S  e  2  is  not  only  homeomorphic  but  also  dif ferentiably  equivalent,  i.e., 
di f  feomorphi c ,  to  S2 .  In  other  words,  topological  transformations  of  S  are 
not  only  continuous  but  smoothly  differentiable.  Notj  that  although 
dif feomorphi sms  are  homeomorphisms,  i.e.,  differentiability  implies 
continuity,  the  converse  is  not  true. 

We  may  summarize  these  considerations  in  the  following. 
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Theorem 


All  the  closed  equipotential  surfaces  2  of  an  idealized  sphere¬ 
shaped  uniformly  rotating  Earth,  i.e.,  those  which  satisfy  the 
Pizzetti  inequality  (3),  are  dif feomorphic  to  a  2-sphere  $2- 
Our  argument  has  been  based  on  translating  the  following  five  physical 
requirements: 

(i)  tne  space  near  the  Earth  is  Euclidean  and  3-dimensional; 

(ii)  each  S  €  2  consists  of  one  piece; 

(iii)  one  can  measure  distances  between  different  points  on  each  S  €  2; 

(iv)  each  S  €  2  admits  am  upward  pointing  direction; 

(v)  neighboring  plumb  lines  on  each  S  6  2  do  not  intersect; 
into  precise  topological  and  geometric  requirements. 

Finally,  the  Gauss-Bonnet  theorem  for  S  yields 

fj~j)  KdA  =  2nx(S)  , 

where  dA  is  the  area  element  on  S  and  *(S)  is  the  Euler-Potncare 
characteristic  on  S.  For  S  topologically  equivalent  to  S2,  x(S)  *  2.  The 
significance  of  this  quantity  is  given  by  the 
Poincare-Hopf  Theorem 

If  v  is  a  smooth  vector  field  on  a  compact  orientable  surface 

S,  and  v  has  only  a  finite  number  of  zeros,  then  the  sum  of  the 

indices  of  singularity  j  of  v  is  equal  to  \(S). 

The  result  \(S)  =  2j  =  2  then  explains  the  inevitability  of  the 
singularities  of  a  smooth  tangent  vector  on  S,  as  mentioned  in  Section  3. 
Typically  one  has  two  singularities  on  S,  viz.,  a  pair  of  sources,  sinks,  or 
centers  with  each  singular  point  having  a  singularity  index  j  =  +1.  The 
diversity  of  possible  singilarities  involving  saddle  points,  dipoles,  etc., 
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but  with  the  sum  of  the  indices  being  constrained  by  \(S)  -  2,  also  indicates 
the  various  choices  of  coordinate  systems  having  ’simpler'  or  more  'severe' 
types  of  singularities. 

The  intent  of  our  discussion  has  been  to  show  that  the  mathematical 
foundations  of  the  Marussi-Hotine  approach  to  geodesy  are  both  rich  and 
intricate.  A  failure  to  succeed  in  resolving  a  particular  question  may  not  be 
attributable  merely  to  a  lack  of  ingenuity  or  initiative  —  a  solution  need 
not  exist!  Lord  Butler  once  said  that  "politics  is  the  art  of  the  possible." 
We  would  suggest  that  a  similar  remark  is  apropos  in  differential  geodesy, 
where  topology  is  the  art  of  the  possible. 
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Figure  4.  The  singularities  of  the  tangent  vectors  on  the  2-schere.  In  the 
left  figure  the  latitudinal  tangents  have  singular  points  at  the  North  and 
South  poles,  each  having  an  index  j  -  +1  corresponding  to  a  "source"  and 
"sink"  respectively  (using  the  language  of  hydrodynamics).  In  the  right 
figure  the  longitudinal  tangents  are  illustrated.  Notice  that  as  the 
geographical  parallels  are  moved  northward  or  southward,  one  obtains 
singularities  again  having  index  j  *»  +1  corresponding  to  "vertices"  at  the 
North  and  South  poles.  The  values  of  the  indices  are  not  direction  dependent! 
This  figure  is  taken  from  Intuitive  Topology  (in  Russian)  by  V.G.  BOLTYANSKII 
and  V.A.  EFREMOVICH,  page  68,  (Nauka,  1982). 
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Figure  5.  Examples  of  singular  points  of  vector  fields,  and  their  indices. 
The  indices  of  these  singularities  are  j  =  +1  in  (a),  (b),  (c);  j  -  -1  in 
(d)f  j  “  -2  in  (e)  and  j  «  +2  in  (f).  This  figure  is  taken  from  H.  HOPF: 
Differential  Geometry  in  the  Large,  Lecture  Notes  in  Mathematics  1000,  page  12 
(Springer-Verlag,  1983). 
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Figure  6.  Examples  of  singular  points  of  vector  fields,  and  their  indices. 
This  shows  that  the  curves,  i.e.  the  trajectories  of  the  vector  fields,  need 
not  be  as  regular-shaped  as  illustrated  in  Figure  5.  This  figure  appears  in 
Topology  from  the  Differential  Viewpoint  by  J.W.  MILNOR,  page  33  (The 
University  Press  of  Virginia,  1965). 


Figure  7.  An  example  of  a  number  of  singular  points  occurring  in  a  single 
chart,  say  on  a  surface  of  irregular  shape.  This  figure  occurs  on  page  73  of 
the  Russian  textbook  cited  in  the  caption  of  Figure  4. 


Appendix 


In  this  Appendix  we  will  illustrate  some  of  the  ideas  in  Sections  4  and  5 
for  the  case  of  a  2-sphere  S2.  Ncno  of  the  results  are  new,  but  the  approach 
may  be  useful  in  showing  how  the  calculation  are  done  when  topological  notions 
are  taken  into  account. 

For  the  2-sphere  S2  the  implicit  equation  is  given  by 
(A.  1 )  FU1,  x2,  x3)  =  (x1)2  +  (x2)2  +  (x3)2  -1  =  0 


where 


denote  Cartesian  coordinates  in  Euclidean  3-space  .  Let 
(A. 2)  Ajk  =  Jl-(xj)2-(xk)2 

denote  the  solution  of  (A.l)  for  one  of  the  coordinates,  say  x1,  in  terms  of 

"i  k 

the  other  two,  viz.  xJ  and  x  with  i,  j  and  k  being  different.  Then 
upon  taking  into  account  the  choice  of  signs  on  the  radical  in  (A.2),  and  the 
fact  one  can  solve  for  the  A^k  in  three  different  ways,  we  obtain  the 
following  six  charts  on  S2: 

(Ulf  x)  =  (x1,  x2,  +A12) 

(U2,  x)  =  (x1,  x2,  -A12* 

(A. 3)  (U3,  x)  =  (x1,  +A13,  x3) 

(Ur,  x)  =  (x1,  -A13,  x3) 

(U^,  x)  =  ( +A 23 r  x2,  x3) 

(U6,  x)  =  (-A23,  x2,  x3)  . 

The  inclusion  of  x  in  the  above  equations  is  necessary  to  indicate  the 
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coordinatization  used  on  the  neighborhoods.  Note  that  the  union  of  the 
coordinate  neighborhoods  U^  and  Uj»  denoted  by  U  covers  S2  except 


for  the  equator 


/  M  ,  2.2  , 

(x  )  +  (x  )  =  1, 


and  it  is  easy  to  see  that  by  forming  the  sixfold  union  U  U  •••  U 
we  obtain  a  covering  of  all  of  $2-  The  six  coordinate  neighborhoods  are 
illustrated  in  Figure  3. 

Using  Hotine's  notation  the  Gaussian  parametrization  on  the  2-sphere  S2 
is  given  by 

(A.  4)  u"  =  (w,  <f>) 


where  <n  is  the  longitude,  and  <p  is  the  latitude  as  illustrated  in  Figure 
4.  We  take  these  to  range  over  the  open  intervals: 


(A. 5) 


0  <  (d  <  2n ,  -tt/2  <  <p  <  tt/2 


and  hence  we  have  excluded  the  meridian  w  «  0  from  our  parametrization. 
Note  that  this  meridian  also  contains  the  North  pole 


and  the  South  pole 


N  :  (<i>  -  0,  4>  -  rr/2) 


S  :  ((.>  =  0,  <t>  =  -tt/2). 


If  we  changed  the  parametrization  (A. 4)  to 


u  m  (<■',  <P ) 


where  c.i  remains  the  longitude,  and  <P  is  the  colatitude,  i.e. 


(A.7 ) 


<t>  =  tt/2  - 


(A. 8) 


0  <  (.>  <  2tT,  0  <  0  <  TT  , 


(A. 9) 


=  -1  . 
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This  sign  is  opposite  that  employed  in  (6),  however  all  that  was  required  is 
that  the  functional  determinant  have  a  fixed  sign. 

r  12  3 

In  terms  of  the  Cartesian  coordinates  x  =  x  =  (x  ,  x  ,  x  )  and  the 
parametrization  (A. 4)  and  (A. 5)  we  have  the  following  parametrization: 


(A. 10) 


1 

x~  =  cos  <•>  cos  <t>  , 
2 

x  =  sin  (.)  cos  <(>  , 
x3  =  sin  <t>  , 


which  by  (A. 5)  omits  the  North  and  South  poles  N  and  S  respectively,  i.e. 
(0,  0,  1)  for  <t>  -  tr/2  and  (0,  0,  -1)  for  <f>  »  -rr/2. 

By  direct  calculation  we  have 


1  2, 

d(x  ,x  ) 

d  (<■>,$) 


=  cos  <t>  sin  <p 


(A. 11) 


2  3 

d(x,x) 

d(w  ,<#>) 


2 

=  cos  <■>  cos 


<p 


t 


3  1. 

d(x  ,x  ) 


sin  w  cos  <P  ; 


d{u,<p) 

and  the  requirements  of  the  Inverse  Function  Theorem  are  satisfied.  Note  that 
at  the  excluded  poles  N  and  S,  all  three  of  these  functional  determinants 
are  identically  zero. 

Thus,  the  port  la l  covering  of  S2  by  coordinate  neighborhoods  (U,  x) 


consists  of  a  pair  of  right  and  left  hemispheres  which  omits  the  meridian  of 
reference  w  =  0. 

Corresponding  to  the  parametrization  (A. 10)  the  tangent  vectors  to  the 
longitude  and  latitude  are  respectively  given  by 
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(A. 12) 


dxC/du  -  (-sin  w  cos  cos  w  cos  $ ,  0) 
dx  /d <t>  =  (-cos  w  sin  4> ,  -sin  u  sin  <#> ,  cos  <t>) . 


We  denote  these  by  x^,  x  respectively,  and  note  that  on  the  equator 


we  have 


x  =  (-sin  w,  cos  (.1,  0) , 


X.  =  (0,  0,  1); 

~9 


however  at  the  excluded  poles  N  and  S 


x  =(0,0,0), 


x  =  (+  cos  c»,  +  sin  ci,  0). 


This  suggests  that  —  at  least  in  this  parametrization  —  both  N  and  S  are 
singular  points  since  the  longitudinal  tangent  vector  x  vanishes  at  these 


points. 


By  using  (10)  we  obtain  the  following  expression  for  the  first 


fundamental  form  (9) 
( A. 13 ) 


2  2  2  2 
ds^  =  cos>  <L>  +  d<T, 


cos  <P  0 


(A. 14] 


(A. 15) 


a  =  cos  <£ 


Note  that  a  >  0,  and  a  =  0  only  at  the  exclude  poles  N  and  S. 
The  covariant  components  of  the  unit  normal  n  are  defined  by 


(A.16: 
and  hence 


1  m  n  nr/3 

II  =  T  t  x  x„  f 
r  2  rmn  a  p 


(A. 17) 


l!  =  (cos  ci  cos  </>,  sin  ci  cos  <t>,  sin  </•). 
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V 


1 


On  the  equator  (w,  <p  =  0)  we  have 

i>  =  (cos  (.),  sin  <<>,  0) 

and  curiously  at  the  excluded  poles  N  and  S 

p  =  (0,  0,  il). 

This  is  surprising  since  as  noted  previously  the  longitudinal  tangent  vanishes 
at  both  N  and  S.  The  result  occurs  since  in  (A.  16)  the  Levi-Civita 
dualizor  ea ^  involves  a  factor  of  l/Va  which  cancels  with  a  common  factor 
of  cos  <#>  in  the  vector  product  of  x  and  x, . 

~<.i  _<f> 

If  we  regard  (A. 10)  as  being  the  imbedding  equations  of  in  E^,  then 

the  range  (A. 5)  of  values  of  m,  <p  omits  not  only  the  reference  meridian 
w  -  0,  but  the  poles  N  and  S.  Hence  (A. 10)  cannot  be  taken  as  defining  all 
of  S2  in  E^.  The  expressions  (A.11)-(A.17)  show  no  analytical  difficulties 
along  the  reference  meridian  except  at  the  excluded  poles  <p  =  ±  v/2. 

of  the  second  fundamental  form 


0 
1 

Note  that  again  b  >  0,  and  b  =  0  only  at  the  excluded  poles  N  and  S.  By 
taking  the  usual  definitions  of  the  Gaussian  (total)  and  Germain  (mean) 
curvature  we  find  that 

K  =  1,  H  =  -1. 

By  direct  calculation  using  the  Gauss-Bonnet  formula,  since  the  area 
element  of  dA  on  is 


Likewise,  computing  the  coefficients  b 


by 

(A. 18) 
we  have 


(A. 19) 


op 


,  K  r  s 

b  n  =  -<5  x  vn 
ap  rs  a  p 


2, 

-  cos  <t> 


ctp 


and 


(A. 20) 


b  =  cos  </> . 
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dA  =  cos  <p  dw  d<£, 


Hence,  x(S 2 )  ■  2  as  claimed.  This  is  a  correct  result,  however  the 

calculation  should  be  done  more  carefully.  We  consider  a  partial  covering  of 

S2  consisting  of  the  Northern  hemisphere 

1  2  2  2  3  2  3 

UN  :  (xV  +  (xV  +  (x  r  =  1,  (x  >  0); 

the  Southern  hemisphere- 

1  2  2  2  3  2  3 

Us  :  (xV  +  (xV  +  (xV  =  1,  (x  <  0); 

and  denote  the  equator  by 

12  2  2 

CQ  :  (xV  +  (x  )  -  1. 

Then  one  calculates 


(A. 21) 


K  dA  + 

JJUN 


where  a  is  the  geodesic  curvature  of  Cq,  ds  is  the  element  of  arc  length 
on  Cq,  and  the  line  integral  is  taken  in  a  sense  such  that  each  of  the 
hemispheres  remains  on  the  left-hand  side  of  CQ.  In  terms  of  the  Gaussian 
parametrization  we  have 


0  <  w  <  2tt  -2tt  <  w  <  0 

(A. 22)  U  :  ;  U  :  " 

0  <  0  <  tt/2  s  -tt/2  <  <p  <  0 


Then  each  of  the  surface  integrals  is  equal  to  2rr,  and  when  the  two 
expressions  in  (A. 21)  are  added,  the  line  integrals  cancel.  Actually  since 
CQ  is  a  great  circle,  then  a  =  0  and  the  integrands  of  each  of  the  line 
integrals  is  identically  zero.  This  can  be  verified  by  using  the  classical 
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expression  for  a: 


(A. 25 ) 
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